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Abstract 

In this paper two common collocation approaches based on radial basis functions have been 
considered; one be computed through the integration process (IRBF) and one be computed 
through the differentiation process (DRBF). We investigated the two approaches on natural 
convection heat transfer equations embedded in porous medium which are of great impor- 
tance in the design of canisters for nuclear wastes disposal. Numerical results show that the 
IRBF be performed much better than the common DRBF, and show good accuracy and 
high rate of convergence of IRBF process. 
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1. Introduction 

Natural convective heat transfer in porous media has received considerable attention 
during the past few decades. This interest can be attributed due to its wide range of 
applications in ceramic processing, nuclear reactor cooling system, crude oil drilling, chemical 
reactor design, ground water pollution and filtration processes. External natural convection 
in aporous medium adjacent to heated bodies was analyzed by Nield and Bejan Merkin 
[2I, [if, Minkowycz and Cheng 0, [sf. Pop and Cheng ja, 1^, Ingham and Pop [s*]. In all of 
these analysis, it was assumed that boundary layer approximations are applicable and the 
coupled set of governing equations were solved by numerical methods. 

In this paper, the same approximations are applied to the problem of natural convection 
about an inverted heated cone embedded in a porous medium of infinite extent. No similarity 
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solution exists for the truncated cone, but for the case of full cone, if the prescribed wall 
temperature or surface heat flux is a power function of distance from the vertex of the 
inverted cone similarity solutions exist [H, loj , a great deal of information is available on heat 
and fluid flow about such cones as reviewed by Refs. j^, [lo| . 

Bejan and Khair used Darcy's law to study the vertical natural convective flows 
driven by temperature and concentration gradients. Nakayama and Hossain 12| applied 
the integral method to obtain the heat and mass transfer by free convection from a vertical 



surface with constant wall temperature and concentration. Yih [l3j examined the coupled 
heat and mass transfer by free convection over a truncated cone in porous media for variable 
wall temperature and concentration or variable heat and mass fluxes and IJ] applied the 
uniform transpiration effect on coupled heat and mass transfer in mixed convection about 
inclined surfaces in porous media for the entire regime. Cheng 15|] used an integral approach 
to study the heat and mass transfer by natural convection from truncated cones in porous 
media with variable wall temperature and studied the Soret and Dufour effects on the 
boundary layer flow due to natural convection heat and mass transfer over a vertical cone in 
a porous medium saturated with Newtonian fluids with constant wall temperature. Natural 
convective mass transfer from upward-pointing vertical cones, embedded in saturated porous 
media, has been studied using the limiting diffusion \v^. The natural convection along 
an isothermal wavy cone embedded in a fluid-saturated porous medium are presented in 
18l [lOj]. Lai and Kulacki 20| studied the natural convection boundary layer flow along a 
vertical surface with constant heat and mass flux including the effect of wall injection. In 



2l| fluid flow and heat transfer of vertical full cone embedded in porous media have been 



solved by Homotopy analysis method j22|, |23 



Mathematical modeling of many problems in science and engineering leads to ordinary 
differential equations (ODEs) [24-28|. The methods based on radial basis functions (RBF) 
which are part of an emerging field of mathematics are famous ways to solve these kinds 
of problems . First studied by Roland Hardy, an Iowa State geodesist, in 1968, these 
methods allow for scattered data to be easily used in computations [29|. The concept of 
solving DEs by using RBFs was first introduced by Kansa |30|. Since then, it has received 
a great deal of attention from researchers. And consequently, many further interesting 
developments and apphcations have been reported (e.g. Zerroukat et al.|3l|, Mai-Duy and 



Tran-Cong 321. |33|). Essentially, in a typical RBF collocation method, each variable and its 
derivatives are all expressed as weighted linear combinations of basis functions, where the 
sets of network weights are identical. These closed forms of representations are substituted 
with the governing equations as well as boundary conditions, and the point collocation 
technique is then employed to discretize the system 27(|. If all basis functions in networks 



are available in analytic forms, the RBF collocation methods can be regarded as truly 



meshless methods 34 



There are two basic approaches for obtaining new basis functions 
from RBFs, namely direct approach (DRBF) based on a differential process (Kansa 301 ]) 
and indirect approach (IRBF) based on an integration process (Mai-Duy and Tran-Cong 



29|, 122, l35|). Both approaches were tested on the solution of second order DEs and the 



indirect approach was found to be superior to the direct approach (Mai-Duy and Tran-Cong 
0). 
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In this paper we apply the DRBF and IRBF for solving natural convection of Darcian 
fluid about a vertical full cone embedded in porous media prescribed surface heat flux which 
is third order nonlinear ODE. 



2. Problem formulation 

Consider an inverted cone with semi-angle 7 and take axes in the manner indicated in 
Fig. [^a). The boundary layer develops over the heated frustum x = xq- 

The boundary layer equations for natural convection of Darcian fluid about a cone are 

-M + -M = o, (2.1) 

^ ^ Poc/3Kg cos7(T - T^) 

dT dT _ d'^T 
dx dy dy"^ 

For a thin boundary layer, r is obtained approximately xsin{'y). Suppose that a power law 
of heat flux is prescribed on the frustum. Accordingly, the boundary conditions at infinity 
are: 

u = 0, T = Too, if y — 00 (2.3) 

(2.4) 

and at the wall are 
v = if y = 0. 
If the surface heat flux \d\ is prescribed, q^ is obtained as 
dT 

Qw = -k{-^)y=o = A{x - xo)^, xo<x<oo. 

For the case of a full cone {xq = 0, Fig^b)) a similarity solution exists joj. 

In the case of prescribed surface heat flux the similarity solution for the stream function 
ip and T where 

Idip Idtp 

u = - — , v = — — 2.5 
r ay r ox 



is of the form 

^ = ar{Ra,y/''fiv), (2.6) 
T-T^=^-^{Ra,)-k{v), 

X 



where 

^ poo/3gKcos{-f)q^x'^ ^2 7) 

is the local Rayleigh number for the case of prescribed surface heat flux. The governing 
equations become 

/' = 0, (2.8) 

2 3 
subjected to boundary conditions as: 

/(0) = 0, ^'(0) = -l, 0(oo) = O. (2.9) 

Finally from Eqs. (12. 8 P and (12. 9p we have: 

(ode. r + (^)//"-m(/')^ = o, 

[B.C. /(O) = 0, no) = -1, /'(oo) = 0. 
It is of interest to obtain the value of the local Nusselt number which is defined as [oj: 

Nu, = (2-11) 

k(T^ TqqJ 

From Eqs. (12.111) . (12. 6p and (12. 7p it follows that the local Nusselt number which is interest 
to obtain given by: 

Nu^ = Ral^'^[-e{0)]. (2.12) 
3. RBF Functions 

Let = {a; G M, a; > 0} be the non-negative half-line and let (p : — )■ M be a 

continuous function with 0(0) > 0. A radial basis function on is a function of the form 



|X-X,1|) (3.1) 

where X, Xj G M'^ and ||.|| denotes the Euclidean distance between X, Xj. If one chooses 
X points {Xi}'^Li in M then by custom 

N 

s(X) = ^A,0(||X-X,||); A, gM, (3.2) 



1=1 



is called a radial basis function as well 36 



In order to explain RBF methods briefly, suppose that the one-dimensional input data point 
set or the center set xt in the given domain f2 C M is given. The center point is not necessarily 
structured, that is, it can have an arbitrary distribution. The arbitrary grid structure is one 
of the major differences between the RBF method and other global methods. Such a mesh- 
free grid structure yields high flexibility especially when the domain is irregular. In this 
work the uniform grid is used for RBF approximation. 
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3.1. Properties of RBF 

With a radial function 0(r) and with data values Ui given at the locations x,, for i = 
1, 2, the function 



TV 



s{x) = \i(f)i{x) (3.3) 



i=l 



where r = = ||x — and (pi{x) = (f){\\x — Xi\\), interpolates the data if we choose the 



expansion coefficients Aj in such a way that s(xj) — tij, for j = 1,2,..., AT [37|, |38[. The 
expansion coefficients Aj can therefore be obtained by solving the linear system AA = U, 
where: 

A, =(j){\\xj-x,\\), (3.4) 
A= [Ai,A2,...,A^]^, (3.5) 

U = [Ui,U2, . . . ,UAr]^. (3.6) 

All the infinitely smooth RBF choices listed in Table ([1]) will give coefficient matrices 



A in ( 13. 4 p which are symmetric and nonsingular |39|, i.e. there is a unique interpolant of 
the form (13. 3 p no matter how the distinct data points are scattered in any number of space 
dimensions. In the cases of inverse quadratic, inverse multiquadric and GA the matrix 
A is positive definite and, for multiquadric (MQ), it has one positive eigenvalue and the 



remaining ones are all negative [39 



Interpolation using Conical splines and thin-plate splines (TPSs) can become singular in 



multidimensions [3J]. However, low-degree polynomials can be added to the RBF interpolant 
to guarantee that the interpolation matrix is positive definite (a stronger condition than 
nonsingularity) . For example, for the Conical RBF and the TPS in d dimensions this 
becomes the case if we use as an interpolant s{x) = J21Li(^iPi{^) + "l^iLi ^i4>{\\^ ~ 
together with the constraints J2j=i ^jPii^j) = O5 ^ = 1,2, ...,m. Here Pi{x) denotes a 
basis for polynomials of in R'^ [V^ denotes the space of d-variate polynomials of order 
not exceeding q) and m = {q — 1 + d)\ / {d\{q — j40|. 

3.2. RBF Interpolation 

One dimensional function u{x) to be interpolated or approximated can be represented 
by an RBF as: 

N 

u{x) ^ s{x) = Xi(j)i{x) = <^^{x)A (3.7) 
1=1 

where 

$^(x) = [0i(x), 02(a;), Mx)], (3.8) 
A = [Ai,A2,...,A;v]^, (3.9) 
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X is the input and AjS are the set of coefficients to be determined. By choosing interpolate 
nodes {xj,j = 1, 2, A^) in f2 U dfl, the function u{x) can be approximated in i7 U dQ. 



N 



Uj = ^Xi(j)i{xj), (j = 1,2, A^). 

i=l 

To brief discussion on coefficient matrix we define: 
AA = U 

where 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



Note that 4>i{xj) = 0(||xj — Xj\\) therefore (pi{xj) = (pj{xi) consequently A = A^ . 
The shape parameter c which is shown in Table ([1]) affects both the accuracy of the approx- 
imation and the conditioning of the interpolation matrix ji^]. In general, for a fixed number 
of N, smaller shape parameters produce the more accurate approximations, but also are 
associated with a poorly conditioned A. The condition number also grows with N for fixed 
values of the shape parameter c. Many researchers (e.g. 41, 4^) have attempted to develop 
algorithms for selecting optimal values of the shape parameter. The optimal choice of the 
shape parameter is still an open question. In practice it is most often selected by brute force. 
Recently, Fornberg et. al. 43|] have developed a Contour-Pade algorithm which is capable of 
stably computing the RBF approximation for all c > ji^]. 



U =[U1,U2, ...,MAr]'^, 

A=[<|.^(a;i),<l>^(x2),...,<f^(x^)r 

/ 01 (Xi) 4>2{Xi) ... 4>n{Xi) \ 

0i(a;2) 02(3^2) ••• 07v(a;2) 

y 0l(xAr) (j)2{xN) ■■■ (pNixN) J 



The following theorem about the convergence of RBF interpolation is discussed [45|, |46 



Theorem 3.1. 



1,2, 



.., A^) are N nodes in Q which is convex, let 



h = max min llx — xJU 

xeQ l<i<N 



when (j){r]) < c(l + \r]\) g^^y y,^^^ satisfies j {u{r]))'^ / (j){r])d7] < 00 we have 



u 



where (f){x) is RBF and the constant c depends on the RBF, d is space dimension, I and 
a are nonnegative integer. R can he seen that not only RBF itself hut also its any order 
derivative has a good convergence. 
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3.3. Direct RBF for ODEs (DRBF) 

In the direct method, the closed form DRBF approximating function fl3.7p is first ob- 
tained from a set of training points, and its derivative of any order, e.g. pth order, can then 
be calculated in a straightforward manner by differentiating such a closed form DRBF as 
follows: 

^ - ^M^)) -t^'W-t ^.«!"(-) (3.14) 

i=l i=l 1=1 

where 

Gf\x) = = 0, 1, ...,p. 

Now we aim to apply the DRBF method for solving the ODEs in general form : 

F(x,u,u',...,u^P-^\u^P'^) = 0, a<x<b, , , 

where Cj G {a,b} and u^^\x) = d^u{x)/dx\ F is known function and are known 

constants. By substituting Eq. (13.71) in fl3.15p and using Eq. fl3.14p we have: 

N N N 

F{x, s, s', = Fix, \Gf\x), ^ XiGf\x), ^ X.G?{x)). 

i=l i=l i=l 

Now, to obtain AjS [i = 1, 2, N) we define the residual function: 

Res{x) = F{x, s, s', s^^^). (3.16) 

The set of equations for obtaining the coefficients {Xi}fL^ come from equalizing Eq. (I3.16P 
to zero at N — p interpolate nodes {xj}fSi plus p boundary conditions: 

(Res{xj) = 0, J = l,2,...,iV-p, ^^^^^ 

Since the direct approach is based on a differentiation process, all derivatives obtained here 
are very sensitive to noise arising from the interpolation of DRBFs from a set of discrete 
data points. Any noise here, even at the small level, will be badly magnified with an increase 



in the order of derivative 29 
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3.4- Indirect RBF for ODEs (IRBF) 

In the indirect method, the formulation of the problem starts with the decomposition 
of the highest order derivative under consideration into RBFs. The obtained derivative 
expression is then integrated to yield expressions for lower order derivatives and finally 
for the original function itself. In contrast, the integration process, where each integral 
represents the area under the corresponding curve, is much less sensitive to noise. Based 
on this observation, it is expected that through the integration process, the approximating 
functions are much smoother and therefore have higher approximation power. Also To 
numerically explore tile IRBF methods with shape parameters for which the interpolation 



matrix is too poorly conditioned to use standard methods [47[. Let p be the highest order 
of the derivative under consideration the boundary value ODEs in general form Eq. fl3.15p 
when (3A; s.t. = h) then we can define: 



dxP 

dP-^s{x) 
dxP-^ 



N 



Y,^M^)^ (3-18) 

1=1 

N AT „ TV 

^^Xi(j)i{x)dx = '^^Xi / (f)i{x)dx = '^^Xih'f ^\x) + di, 

i=l i=l i=l 



^ = / E ^^hf{x)dx + d,j^_ + d,j^_ + ... + 



i=l 

N 



„p-2 ry.p-'i 



S[X] 



J (2^ ^ih\ \x) + di + ^2 (p_3)i + - + dp-i)dx 

2^ Xih\ '(X) + d, + '^^ {p_2)\ ^ '^"-^'^ ^ 



where 



\jh\ \x)dx, k = 0,1, ...,p — 2. 

Substituting Eqs. (13.181) in fl3.16p at interpolate nodes {xj}jLi plus p boundary conditions 
the set of coefficients {Xi}fL^ and {dj}^^^ is obtained as follow : 

'Res{xj) = 0, J = 1,2,..., AT, 
sW(ei) = ai+i, i = 1,2, ...,p. 



4. Solving the model 



Consider governing equation of fluid flow and heat transfer of full cone embedded in 
porous medium that is expressed by Eq. (12.101) for prescribed surface heat flux. 

In the flrst step of our analysis, we approximate f{ri) for solving the model by DRBF: 



N 



f{v)^s{v) = J2xMv), (4.1) 



1=1 



and f"'{r]) for solving the model by IRBF: 

r(^)=^^S(r^) = E^^'^^(^)- (4-2) 



1=1 



The general form of problem appear to: 



/(0)=0, r(0) = -l, /'(oo) = 0. 



(4.3) 



To solve this problem we deflne residual function: 

Res{7]) = Fxiv, s, s', s", s'"); for DRBF, (4.4) 

Res{7]) = Fxiv, s, s, s", s'"); for IRBF. (4.5) 

The unknown coefficients {Xi}f^i come from equalizing Res{r]) to zero at N interpolate 
nodes rji from Uniform distribution between and r]^ which we set 9/2 for this problem. 

4.I. Solving the model by DRBF 

In the flrst step of solving, (piirj) is set by inverse multiquadric function which is shown in 
Table ([T]). Now, the residual function is constructed by substituting Eq. (14. ip in Eq. (14. 4p : 

N . TV N 9\ 1 ^ 

1=1 i=l 1=1 1=1 

By using N — 2 interpolate nodes {r]j}^Si plus two boundary conditions of Eq. (14. 3 p 
(/(O) = 0, /"(O) = —1), the set of equations can be solved, consequently the coefficients 
will be obtained: 

'Resir],)=0, j = 1, 2, iV - 2, 

EiiA.0.(O) = O, 

.Ef=iA'<^^'(0) = -l- 
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Take into account (pi{c>o) = 0, for i = 1,2, ...,N the infinity boundary condition (/'(oo) = 0) 
is already satisfied. 



Table (|2]) show the f'{ri) for some A in comparison with solutions of [2l|]. Also f'{r]) for two 
selected A = 1/4 and 3/4 are showed in Table (j3]) in comparison with Runge-Kutta solution 
is obtained by the MATLAB software command ODE45 which is used and applied by the 
authors in ref. j2l|. Absolute errors show that DRBF give us approximate solution with 
a high degree of accuracy with a small A^. The resulting graph of Eq. fl2.10p is shown in 
Figure 

4-2. Solving the model by IRBF 

In the first of solving 0i(?7) is set by multiquadric function which is shown in Table (JT]). 
Now, the residual function is constructed by substituting Eq. (14.21) in Eq. (14. 5 p and using 

Eq. fiarrsj) : 

^ . ^ N ^ ^ ^ N . 

(f)i{r])dr] 



i=l i=l J J J j=i J 



i=l i=l i=l 

,2A + 1, ^ 



-)C£^^h^^^^) + dlV + d 



2 

2) ■ 



i=l 

By using A^ interpolate nodes {r/jjjLi plus three boundary conditions of Eq. (14.31) the set of 
equations can be solved, consequently the coefficients {Xi}fLi and {rfj}^^]^ will be obtained. 
In this method we put r]^ instead of infinity condition: 

(Res{T]j) = 0, J = 1,2,..., AT, 

Ef=iA./i?'(0) + rf3 = 0, 

E|iMf(o) + rfi = -i, 



Table (jl]) shows the /'{r]) for some A in comparison with solutions of [2l|. Also /'(?]) 
for two selected A = 1/4 and 3/4 are showed in Table (E]) in comparison with Runge-Kutta 
solution is obtained by the MATLAB software command ODE45 which is used and applied 
by the authors in ref. jllj. Absolute errors show that IRBF give us approximate solution 
with a high degree of accuracy with a small A^. The resulting graph of Eq. (I2.10p is shown 
in Figure ([3]). A graph in figures (jl]) for A = 2/3 show ||i?es||^ for some A^. Table ([6]) shows 
||-Res||^ for some A^ and A. 

Comparison between DRBF solution in Table ([3]) with A^ = 12 and IRBF solution in Table 
(E]) with A^ = 10 for /'(//) show that the convergence of the IRBF method is faster, because 
of using less numbers of collocation points. 
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5. Conclusion 

In this paper we made a comparison between the two common collocation approaches 
based on radial basis functions namely DRBF and IRBF methods on natural convection 
equation about an inverted heated cone embedded in a porous medium of infinite extent 
which are of great importance in the design of canisters for nuclear wastes disposal. These 
functions are proposed to provide an effective but simple way to improve the convergence 
of the solution by collocation method. The direct approach (DRBF) is based on a differen- 
tiation process, all derivatives are very sensitive to noise arising from the interpolation of 
DRBFs from a set of discrete data points. Any noise, even at the small level, will be bad 
magnified with an increase in the order of derivative. The indirect technique (IRBF) which 
is based on integration process, each integral represents the area under the corresponding 
curve, is much less sensitive to noise. Based on this observation, it is expected that through 
the integration process, the approximating functions are much smoother and therefore have 
higher approximation power. Additionally, through the comparison with other methods 
such as HAM we show that the RBFs methods have good reliability and efficiency. Also 
high convergence rates and good accuracy are obtained with the proposed method using 
relatively low numbers of data points. 
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Figure 2: DRBF approximation of /'(r/) for different values A = 0, 1/4, 1/3, 1/2, 3/4 and 1 
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Figure 3: IRBF approximation of f'{r]) for different values A = 0, 1/4, 1/3, 1/2, 3/4 and 1 
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Figure 4: ||-Res||^ for A = 
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Nomenclature 


A 


prescribed constant 


f 


similarity function for stream function temperature 


9 


acceleration due to gravity parameter 


K 


permeability of the fluid-saturated porous medium 




local Nusselt number 


Qw 


surface heat flux fluid-saturated porous medium 


r 


local radius of the cone fluid 


RcLx 


local Raleigh number 


T 


temperature 


T 


ambient temperature 


U, V 


velocity vector along x,y axis 


x,y 


Cartesian coordinate system 


Xq 


distance of start point of cone from the vertex 


Greek symbols 


a 


thermal diffusivity the fluid-saturated porous medium 




expansion coefficient of the fluid 


V 


independent dimensionless 


e 


similarity function for 


A 


prescribed constants 




viscosity of the fluid 


Poo 


density of the fluid at infinity 




stream function 
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Table 1: Some well-known functions that generate RBFs {r = \\x — Xi\\ = ri), c > 



Name of functions Definition 
Multiquadrics (MQ) + 

Inverse multiquadrics (IMQ) l/{\/r'^ + c^) 
Thin plate (polyharmonic) Splines (TPS) [—l)^^^r'^^log{r) 

Conical splines ^2fe+i 

Gaussian (GA) exp{—cr'^) 

Exponential spline exp{—cr) 
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Table 2: A comparison between solutions obtained by 



21[ and the DRBF method for f'jO) 





Runge-Kutta 






DRBF method 




Other methods 


A 


Solution[21] 


N 


c 


DRBF 


Error with RK 


HAM [21] 





0.94760 


10 


3.46543 


0.94750 


0.0001 


0.94783 


1/4 


0.91130 


12 


3.943 


0.91086 


0.00044 


0.91119 


1/3 


0.90030 


10 


4.9665 


0.90038 


0.0008 


0.90103 


1/2 


0.87980 


10 


5.36 


0.87981 


0.00001 


0.87964 


3/4 


0.85220 


12 


5.23 


0.85227 


0.00007 


0.85242 


1 


0.82760 


10 


5.89 


0.82737 


0.00023 


0.82726 



Table 3: Comparison between DRBF solution and Runge-Kutta solution for /'(??) for A = 1/4 and A = 3/4 

with iV = 12 

A = 1/4, /'(r/) A = 3/4, f (r/) 



DRBF Runge-Kutta Error DRBF Runge-Kutta Error 

rj Solution Solution[21] with RK Solution Solution[21] with RK 









910886 





911295 





000409 





852268 





852193 


0.000075 


0.1 





813122 





813604 





000482 





755678 





755377 


0.000301 


0.2 





720720 





721351 





000631 





666229 





665448 


0.000781 


0.3 





634656 





635531 





000875 





584176 





582985 


0.001191 


0.4 





555601 





556661 





001061 





509628 





508141 


0.001487 


0.5 





483877 





484997 





001120 





442519 





440849 


0.001670 


0.6 





419467 





420587 





001120 





382620 





380907 


0.001713 


0.7 





362196 





363276 





001080 





329539 





327973 


0.001566 


0.8 





311655 





312677 





001022 





282878 





281536 


0.001342 


0.9 





267358 





268264 





000906 





242076 





241013 


0.001043 


1 





228756 





229508 





000752 





206652 





205832 


0.000820 


1.1 





195268 





195878 





000610 





176015 





175434 


0.000581 


1.2 





166342 





166847 





000505 





149653 





149275 


0.000378 


1.3 





141450 





141837 





000387 





127026 





126821 


0.000205 


1.4 





120090 





120362 





000272 





107681 





107596 


0.000085 


1.5 





101811 





102025 





000214 





091167 





091196 


0.000029 
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Table 4: A comparison between solutions obtained by |21[ and the IRBF method for /'(O) 



Runge-Kutta 
A Solution [2JJ 



IRBF method 



Other methods 



N 



IRBF 



Error with RK 



HAMf21 





1/4 
1/3 
1/2 
3/4 
1 



0.94760 
0.91130 
0.90030 
0.87980 
0.85220 
0.82760 



10 
10 
10 
10 
10 
10 



1.860 
2.005 
2.050 
2.150 
2.418 
2.380 



0.94758 
0.91128 
0.90030 
0.87979 
0.85206 
0.82762 



0.00002 
0.00002 
0.00000 
0.00001 
0.00014 
0.00002 



0.94783 
0.91119 
0.90103 
0.87964 
0.85242 
0.82726 



Table 5: Comparison between IRBF solution and Runge-Kutta solution for f'{ri) for X — 1/4 and A = 3/4 



with N 


= 10 














A = 1/4, f'iv) 












IRBF 


Runge-Kutta 


=; 

Error 


LSXDr 


t5 t7 

tiunge-Kutta 


==; 

Error 


V 


Solution 


Solution[21] 


with RK 


Solution 


Solution [21] 


with RK 





0.911278 


0.911295 


0.000017 


0.852059 


0.852193 


0.000134 


0.1 


0.813594 


0.813604 


0.000010 


0.755254 


0.755377 


0.000123 


0.2 


0.721302 


0.721351 


0.000049 


0.665278 


0.665448 


0.000170 


0.3 


0.635394 


0.635531 


0.000137 


0.582736 


0.582985 


0.000249 


0.4 


0.556440 


0.556661 


0.000221 


0.507827 


0.508141 


0.000314 


0.5 


0.484753 


0.484997 


0.000244 


0.440510 


0.440849 


0.000339 


0.6 


0.420334 


0.420587 


0.000253 


0.380599 


0.380907 


0.000308 


0.7 


0.362989 


0.363276 


0.000287 


0.327647 


0.327973 


0.000326 


0.8 


0.312347 


0.312677 


0.000330 


0.281184 


0.281536 


0.000352 


0.9 


0.267934 


0.268264 


0.000330 


0.240696 


0.241013 


0.000317 


1 


0.229197 


0.229508 


0.000311 


0.205482 


0.205832 


0.000350 


1.1 


0.195584 


0.195878 


0.000294 


0.175053 


0.175434 


0.000381 


1.2 


0.166546 


0.166847 


0.000301 


0.148891 


0.149275 


0.000384 


1.3 


0.141548 


0.141837 


0.000289 


0.126408 


0.126821 


0.000413 


1.4 


0.120107 


0.120362 


0.000255 


0.107158 


0.107596 


0.000438 


1.5 


0.101771 


0.102025 


0.000254 


0.090757 


0.091196 


0.000439 


2 


0.043769 


0.043951 


0.000182 


0.038887 


0.039223 


0.000336 


2.5 


0.018310 


0.018546 


0.000236 


0.016225 


0.016574 


0.000349 


3 


0.007391 


0.007610 


0.000219 


0.006433 


0.006832 


0.000399 


3.5 


0.002716 


0.002953 


0.000237 


0.002255 


0.002668 


0.000413 


4 


0.000719 


0.000962 


0.000243 


0.000445 


0.000913 


0.000468 


4.5 


0.000010 


0.000123 


0.000113 


-0.00001 


0.000237 


0.000247 
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Table 6: 


\\Resf for different N and A by IRBF 








A 


N=5 


N=6 


N=8 


N=10 


N=12 




N=15 







1.161037 


0.306257 


0.016229 


0.000108 


0.77e - 


5 


O.lle - 


7 


1/4 


0.851976 


0.278068 


0.018364 


0.000363 


0.52e- 


4 


0.27e - 


6 


1/3 


0.770207 


0.269714 


0.018842 


0.000494 


0.24e - 


4 


0.13e- 


5 


1/2 


0.637068 


0.216095 


0.016813 


0.000587 


0.83e - 


5 


0.47e - 


6 


3/4 


0.469112 


0.230425 


0.020280 


0.001051 


0.15e- 


4 


0.83e - 


6 


1 


0.355474 


0.221336 


0.020463 


0.001077 


0.31e- 


4 


0.39e - 


6 
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